Statistical mechanics and path integrals for a finite number of bosons. 
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Recent investigations show that the statistical mechanics of a finite number of particles in ideal 
harmonic systems predicts different results for the same physical properties, depending on the en- 
semble under consideration. Path integral methods for a finite number of bosons with equidistant 
energy levels give the same answers for the mean energy, the specific heat and the condensation 
temperature etc., irrespective whether their calculation results from the density of states, from the 
partition function or from the generating function. 

We show that this contradiction is due either to the use of approximate relations between quantum 
statistical expressions, or to a misinterpretation of the generating function. 
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Since Bose-condensed vapors became available to experimentation |l| ^] , considerable theoretical interest has been 
raised in the statistical mechanics of an idealization of the vapor, i.e. a system with equally spaced energy levels for 
a finite number of particles. Conventional statistical mechanics for a system in equilibrium at a given temperature 
(expressed by (3 = 1/kT ) states that the probability of the system having energy E n is proportional to exp (—/3E n ) . 
The calculation of the energies E n is a quantum mechanical problem. Equally spaced energy levels are justified by the 
parabolic confinement potential. The quantum statistical theory takes into account the discreteness of the levels and 
the limited number of particles and contrasts therefore with earlier studies of Bose-Einstein condensation, because 
the thermodynamic limit and the quasi-continuity of the levels cannot be used as a justification ||. In the path 
integral approach to quantum statistical theory ||[7| , the calculation of the levels as well as the density of states or 
the partition function form an intrinsic part of the same study and are obtained simultaneously, in contradistinction 
with the approach where a quantum problem gives the energy levels and where subsequently statistical theory is used 
to study the cooperative behavior. 

For identical particles the conditional probability to find a system containing N particles in the neighborhood of the 
configuration r can be calculated by the path integral method We have calculated the partition function Z (f3, N) 
and some static response functions of a Gaussian model for bosons as well as for fermions using this approach. The 
related thermodynamical quantities have been obtained from the free energy that was calculated using the probability 
generating function S (/3, 7). The partition function Z (f3, N) relates to the density of states ft (E, N) as follows: 



Z(0,N) = J e- /3E n{E,N)dE, (1) 
while the generating function is defined by: 

S(/3,7) = EV^(/W (2) 

N 

It should be remarked that the real variables 7 and (3 used to interrelate the three quantities (E, N) , Z (/?, N) and 



S (/?, 7) are independent variables 10 1. The inversion formula to obtain {E, N) from Z (/3, N) is the inverse Laplace 
transform 

0+ioc 

' ' e TE ' 

(3—ioQ 



Sl(E,N) = — / e T "Z{T,N)dT, (3) 
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where (3 is larger than the real part of all poles of Z (r, N) in the complex r plane. For the inversion of 5 (j3, 7) the 
residue theorem can be used: 

Z ^ = hj^ dz - (4) 

It is generally accepted that the probabilities associated with events by statistical mechanics, are consistent with 
those of quantum statistical theory Referring to the normalization factor Q = exp (—/3E n ) of the probability 
in statistical mechanics as the canonical partition function, it is implicitly assumed that Q and Z are the same for 
the same system, and lead consequently to the same predictions. Therefore it should not matter which one of the 
three quantities Q, Z or S is calculated because of their exact interrelationship. Quantum statistical theory would be 
in conflict with statistical mechanics if in statistical mechanics the quantities corresponding with Q, Z and S (usually 
indicated by their ensemble: microcanonical, canonical or grand canonical), would lead to different predictions, what 
they do according to p2 ]. Recently this puzzle even gave rise to the introduction of a new fourth ensemble, coined the 
"Maxwell Demon" ensemble |l^,|l4|, with the aim to take the ensemble conditions for condensation better into account. 
The quantum statistical analogue of this "Maxwell Demon" ensemble is obtained by performing the inversions starting 
from 5 in the appropriate order 

/3+ioo 

T(B, 7 ) = -L / e^S(r, 7 )dr = ^fi(i?,7V) 7 JV , (5) 

a ■ N 
p—100 

expressing also this quantity in a unique way in terms of the density of states. 

The main question that will be addressed in the present letter is: what is the origin of the apparent discrepancies 
between predictions based on statistical mechanics ? We will discuss here two possibilities. The first one is based on 
probability considerations. The second one elaborates upon an old warning of Zip, Uhlenbeck and Kac jllj against 
grand canonical partition functions calculated from S f/3, e^ M ) = Tr [exp (— (3 (H — fiN))] . 



1. A conditional probability density approach. 

First we show how the discrepancies mentioned above are induced by an approximation which is justified in the 
limit of an arbitrarily large number of particles, but not for a finite number of particles. We do this on the basis of 
our approach B that we briefly describe here. Denoting by r the configuration (ri,r2, . . . , rjv) of the N particles, the 
potential energy function V (r) of our model is given by 



JV 



i=i 



n N 

moj \ 



j<k 



(6) 



Because the model is Gaussian the propagator for the particles (considered distinguishable) is known. Using path 
integral methods the projection on the symmetric boson states (antisymmetric for fermions) could be performed, 
leading to an expression for the partition function Zg ((3, N) for bosons. The center-of-mass contribution factorizes 
out, leaving 3 (JV — 1) internal degrees of freedom with frequency w = V^ 2 — Nu> 2 . Introducing b — e~^ nw ', the 
generating function for these internal degrees of freedom can be written as 



3b (7, 0) = ex P 



t 



e=i 



1 — bt 



(7) 



and is the result of a functional integration. Another but more familiar representation of the same function is: 

1(^+1)^+2) 



s*(7,/j)=n(i 



(8) 



The same function can be obtained directly using combinatorial analysis within the assumptions of statistical me- 
chanics. Jl^-pOj] . In our case (7, f3) is only formally the grand canonical partition function of a set of identical 
particles in a parabolic well. To become really a grand canonical partition function the substitution 7 — > exp(^/3) has 
to be made and a calculation procedure for [i has to be given. If N is sufficiently large, stationary phase or steepest 
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descent methods can be used to approximate the inversion formulas ]ic| ] and give a meaning to the chemical potential 

As an illustration, let us consider the partition function Zb {(3, N) based on the contour integral (Q) for the specific 
example of harmonically interacting bosons in an harmonic confinement potential, i.e. the actual calculations are 
performed using the explicit form (Q). The generating function S# (z, (3) makes the direct numerical evaluation of (^) 
unfeasible. The interested reader is referred to for a discussion. However, considering a circular contour around 
the origin with radius u, the substitution z — ue transforms the contour integral into 

Z B (/?, N) = A- E ( U6% yh - im d9 = A r ^(•""^-"He-^dfl. (9) 



The extremum of [in 2 (ue %6 ,(3) — Nlnu] on the real axis satisfies the condition 

N = u-^- In 3 {u, (3) , (10) 
du 

which is precisely the expression for the expected number of particles N in the "grand canonical ensemble" if u is 
interpreted as u = e^ M . Factorizing out this steepest descent contribution and extracting the real part, one obtains 
with little effort 

Z B {f3, N) = 4 0) 09, N) r * (6) d6, (11) 



ii 



zf IJ>,N) =^2l, (12) 

(13) 

7T ^B(U,f3) 



The corresponding free energy thus becomes the sum of two contributions 

,(0) 1 



F B {(3,N) =F B ">(f3,N)- -In [ \ *(9)d9), (1.4) 



Fg>((3,N)=-~ln^y^. (15) 

Fg" 1 (/3, N) is the result which one would obtain in the "grand canonical treatment" . The remainder i In (J^ 'J' (0) dO) 

can be obtained by integration and this correction is crucial for a finite number of particles. The integrand ^ (9) is 
shown in Fig. 1 for N = 10 as a function of 9 for various temperatures T, expressed in units of the condensation 

1 /3 

temperature T c = (N/((3)) ' fiw/k. Below the condensation temperature the oscillations in ^ (9) are distributed 
more or less uniformly (for bosons) over the integration interval [0, n] and above the condensation temperature they 
are strongly damped. The difference between Fb {13, N) and the zero-order approximation F^ (/3, N) decreases with 
an increasing number of particles. The results for the free energy are shown in Fig. 2-4 for N = 1, 10 and 100, where 
/ {(3, N) = Fb {(3, N) /Nhw is plotted versus T/T c . For comparison, the zero-order approximation F B °^ (/3, N) /Nfiw 
is also plotted. The numerical results for F B (/3, N) are to within 6 digits in agreement with those obtained earlier Q 
from a recurrence relation for the partition function. 

If the partition function obtained using our path integral method relates to the conditional probability function 
for the energy density given the number of particles, we may summarize that the predictions based on the canonical 
partition function in statistical mechanics are not in conflict with quantum statistical theory, provided the density 
of states (microcanonical) is correctly derived by the inversion formulas. This concludes our discussion based on the 
conditional probability interpretation of the density of states. 



2. The joint probability density approach 

It is a common error to interpret a conditional probability as a joint probability function. In order to analyse the 
consequences of this possibility, let us assume that the function (7, (3) is the generator of such a joint probability 
function for the energy E and the number particles N , i.e. we give ourself a marginal distribution for the number of 
particles consistent with that generating function. The theory of continuous-time Markov processes for a queue with 
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infinitely many servers |2lj | might illustrate this argument. In this stochastic model of a facility, customers are entering 
according to a Poisson distribution and are subsequently served with an exponentially distributed service time. In the 
case that both distributions (of the customers and of the servers) have the same parameter, the probability pjv(i) that 
N customers are waiting to be served at time t can be calculated from the generator A of the process. This generator 
explicitly takes the transitions to one more or to one less customer into account. Denoting pjv(i) = (N \ e tA \ iV) and 
g (z, t) — X)jv=o zN PN(t), the analogy between the probability generating function of the stochastic process and the 
generating function for the partition function becomes obvious. 

We will discuss the particle fluctuations of the "grand canonical ensemble" in this interpretation and elucidate the 
remark of ref. [p"lf earlier cited on that ensemble. In order to do so it is instructive to introduce an alternative to invert 
a power series. This inversion will allow to make a connection between the generating function in this interpretation 
and coherent states. Using 
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-W 2 z N z L d 2 z = 5 N>L , (16) 



where z denotes the complex conjugate of z, one obtains the following expression for the partition function Z (/3, N) 
from S (z, (3): 

Z((3,N) = 4r- / e-W 2 z N Z(z,(3)d 2 z. (17) 



Nl 7T _ 

Filling out the power series (|J) and summing over N one finds: 

oo oo oo N y L p -i\z\ 2 

N=Q J N=0L=0 V V 

where use has been made of the orthogonality relation (|l^) to replace a denominator TV! by V N\L\. 

Introducing formally a Hamiltonian TL that allows for transitions between systems with different numbers of parti- 
cles, the following form 

oo . 

J2Z(P,N) = - / (z\e-P H \z)d 2 z (19) 

N=0 n J 

helps to recognize the coherent state representation for the normalization of the probability for the event that "the 
system contains N particles given the inverse temperature f3" . The state | z) can be considered as a coherent matter 
state built up from Bose systems containing a different number of particles: 



E 



p 2 l z l 



v AP 

N v 



\N) . (20) 



The derivation of the actual form of the Hamiltonian TL needed to study other than equilibrium properties is beyond 
the scope of this letter. In order that the equilibrium properties predicted by TL coincide with those derived from 

the generating function, TL should satisfy the condition J^l (N \H\ L) yj^r — 0, which expresses the conservation 
of probability [p2|]25| l, for a generator of a stochastic process. When the equilibrium properties of a model with 
fluctuations in the number of particles are considered, information on this model can be extracted from the generating 
function 5 (z, (3) by making the following identification 

(z\e- l3n \z)=e-^-^S(z,/3). (21) 

The analogy with the queue suggests that if the generating function derives from a density of states that is interpreted 
as a joint probability density, it leads to a marginal distribution of the fluctuating number of particles. Also the 
assumed stochastic behavior requires a specific form for the transitions between a well with A^ bosons and one with 
N ± 1 bosons. These "a priori" properties of that model are certainly not satisfied in the experimental set-up [fi]-j3|]2~i|l 
where the number of particles are monitored. Therefore, interpreting S (z, (3) given by (0) as the generator of a joint 
probability density does in our opinion not allow to draw conclusions on the Bose condensed systems, discussed in 
the introduction. Changing the order of inversion to obtain the "Maxwell demon" distribution without changing the 
interpretation will not alter our objections. This means that the study of a Bose condensed system with respect to a 
fluctuating number of particles requires a quantum statistical theory incorporating explicitly the transitions between 
states containing a different number of bosons, what is not the case in the quantum statistical theory leading to (m). 
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3. Conclusion 



In general we conclude that for a given number of particles the concepts put forward in Q and worked out partially 
in and |)| for systems of identical particles with equally spaced energy levels, allow to trace back differences linked 
with the statistical ensembles to different approximations in the inversion of the generating function. These differences 
have been quantified for the free energy and are identified with oscillatory behavior in our integrand '5 (6) . At least for 
indistinguishable (identical) particles, these oscillations become negligible in the large N limit. It should be stressed 
that the large N behavior is a direct consequence of the projection on the symmetric (antisymmetric) representation of 
the permutation group ensuring indistinguishability of the identical particles. Without the projection on the symmetric 
irreducible representation, this large N behavior for the oscillations is absent. In that case the relative error made 
by replacing the integral by its steepest descent approximation does not become negligible. This observation raises 
a new challenging question in the statistical description of identical particles: is the projection on the symmetric 
(antisymmetric) representation necessary in all circumstances? 
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Figure captions 

Fig. 1: Integrand \& (9) [see for the remainder term of the partition function, calculated for 10 bosons at various 
temperatures. 

Fig. 2: Scaled free energy / (P,N) = Fb ((3,N) /Nhw as a function of T/T c for N = 1. For comparison, the zero- 
order steepest descent approximation F B °^ (P,N) /Nhw is also plotted (dashed line). 

Fig. 3: Same as Fig. 2, for N — 10 bosons. 

Fig. 4: Same as Fig. 3, for N — 100 bosons. 
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